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Abstract. We cross-correlate the 100 um Improved Reprocessing of the IRAS Survey (IRIS) 
map and galaxy clusters at 0.1 < z < 0.3 in the maxBCG catalogue taken from the Sloan 
Digital Sky Survey, measuring an angular cross-power spectrum over multipole moments 
150 < £ < 3000 at a total significance of over 40o\ The cross-spectrum, which arises from 
the spatial correlation between unresolved dusty galaxies that make up the cosmic infrared 
background (CIB) in the IRIS map and the galaxy clusters, is well-fit by a single power 
law with an index of —1.28 ± 0.12, similar to the clustering of unresolved galaxies from 
cross-correlating far-infrared and submillimetre maps at longer wavelengths. Using a recent, 
phenomenological model for the spectral and clustering properties of the IRIS galaxies, we 
constrain the large-scale bias of the maxBCG clusters to be 2.6 ±1.4, consistent with existing 
analyses of the real-space cluster correlation function. The success of our method suggests 
that future CIB-optical cross-correlations using Planck and Herschel data will significantly 
improve our understanding of the clustering and redshift distribution of the faint CIB sources. 
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1 Introduction 

The cosmic infrared background (CIB), which peaks at ~ 200 |xm, largely arises from thermal 
re-emission of UV and optical starlight by dust grains, and comprises nearly half of all 
extragalactic radiation [1-3]. Because ground-based observations suffer from atmospheric 
absorption and balloon/space-based telescopes have had limited resolution, characterising 
the far-infrared and submillimetre (submm) CIB using individual sources — predominantly 
dusty, star-forming galaxies (DSFGs) at < z < 4 — has only been possible fairly recently by 
using statistical techniques, such as stacking (e.g., [4-8]). 

Given the difficulty in resolving CIB sources at long wavelengths, studying fluctuations 
in the CIB (e.g., [9-11]) has emerged as a complementary method to constrain its properties. 
Clustering of unresolved CIB sources was first detected at 160 [im in Spitzer data [12, 13], 
and has since been measured by the Balloon-borne Large Aperture Submillimeter Telescope 
(BLAST; [14, 15]), Herschel's Spectral and Photometric Imaging Receiver (SPIRE; e.g., [16, 
17]) and Planck's High Frequency Instrument [18]. The South Pole Telescope [19] detected 
clustering of dusty galaxies at millimetre wavelengths; shortly thereafter, Ref. [15] measured 
the clustering of the CIB by cross-correlating submm BLAST maps with millimetre Atacama 
Cosmology Telescope maps. 

The clustering of unresolved CIB sources has been used to infer their bias relative to 
the underlying matter fluctuations and to characterise the dark matter halos with which 
they are associated (e.g., [14, 16, 18, 20, 21]). However, extracting clustering properties from 
angular power spectra of CIB anisotropics is challenging because of Galactic dust emission on 
large angular scales and shot-noise from auto-correlation of unresolved sources on small an- 
gular scales. The redshift-distribution of the unresolved CIB sources also remains uncertain, 
particularly at wavelengths longward of 500 |Jm. 

One way to reduce or bypass these issues is to cross-correlate maps containing unresolved 
CIB sources with a catalogue of extragalactic sources with known redshifts. We would not 
expect Galactic emission to be correlated with such a catalogue, so while Galactic dust exists 
as effective noise in the CIB map, it would not add to the cross-correlation signal. A shot- 
noise term on small scales would indicate the extent to which CIB emission is associated with 
the catalogue members, which could be of interest in its own right. 
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In this paper, we present such an analysis by cross-correlating the lOOum Improved 
Reprocessing of the IRAS Survey (IRIS) map [22] with the positions of optically-selected 
galaxy clusters in the maxBCG catalogue [23]. Clustering of the DSFGs in IRIS map has 
been studied before. Ref. [15] marginally detected a weak clustering term in the IRIS power 
spectrum of a ~ 9 deg 2 field with minimum emission from the Galactic cirrus near the South 
Ecliptic Pole. Later, Ref. [24] measured the CIB anisotropics in IRIS maps that were cleaned 
of Galactic cirrus using cross-correlation with HI maps. Our measurement is not sensitive 
to assumptions about Galactic contamination, which acts only as additional noise, as men- 
tioned above, and is not dependent on any external tracer of the Galactic emission; the only 
significant contribution to the cross-spectrum signal comes from the spatial correlation of the 
IRIS DSFGs and the maxBCG clusters. 

The spatial clustering of DSFGs and galaxy clusters is also relevant for ongoing mea- 
surements of the kinematic Sunyaev-Zel'dovich (kSZ; [25]) effect, i.e., the Doppler shift of 
cosmic microwave background (CMB) photons by electrons possessing some bulk motion 
in the presence of a density or ionisation inhomogeneity. Measuring the kSZ angular power 
spectrum is of great interest for constraining an extended, 'patchy' reionisation scenario (e.g., 
[26-28]), but current constraints are degraded because of uncertainty in the degree of corre- 
lation between the CIB sources and the thermal Sunyaev Zel'dovich effect (tSZ; [29]) from 
the scattering of CMB photons off energetic electrons in the potential wells of galaxy groups 
and clusters (e.g., [30-33]). 

This paper is organised as follows. In Section 2, we describe our measurement of the 
IRISxmaxBCG cross-spectrum and our uncertainty estimates; then, in Section 3, we describe 
and discuss a simple model that we fit to the measured spectrum. A discussion and conclusion 
follow in Sections 4 and 5. 

2 Data Reduction 

2.1 Data Sets, Maps and Masks 

Here, we describe the two data sets used in this paper: the IRIS 100 urn map and the 
maxBCG catalogue. 

The IRIS maps are a reprocessed version of the Infrared Astronomical Satellite (IRAS) 
survey [34] which mapped 98% of the sky at 12, 25, 60 and 100 urn. Though IRAS was 
originally designed for discovering and characterising compact sources, it was increasingly 
found useful for studying diffuse, infrared emission. To aid in the latter, particularly with 
the goal of having suitable maps for foreground cleaning of the upcoming Planck satellite, 
IRIS reprocessed the IRAS plates, significantly improving large-scale calibration, destriping 
and zodiacal light removal. In this paper, we use their publicly-available full-sky, 100 urn, 
'no hole' map with HEALPix [35] resolution of NSIDE = 1024. 1 It has a noise level of 
(0.06±0.02) MJy sr" 1 and an effective resolution of (4.3±0.2) arcmin, with the latter obtained 
by assuming that the effective point spread function was a Gaussian. The overall calibration 
uncertainty is estimated at 13.5% and is limited by the uncertainty of the Diffuse Infrared 
Background Experiment (DIRBE) maps which were used for calibration. More details on 
the IRIS maps are in Ref. [22]. 

The maxBCG catalogue of galaxy clusters was compiled by running the maxBCG red- 
sequence selection algorithm [23] on the Sloan Digital Sky Survey (SDSS) photometric data 

: http: //www . cita.utoronto . ca/~mamd/IRIS/Ir isDownload.html 
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Figure 1. Comparison of the redshift distributions of the maxBCG cluster catalogue and the un- 
resolved IRIS galaxies. The solid red histogram shows the distribution of maxBCG clusters, dN/dz 
(with dz = 0.01), where redshifts are photometric and have errors < 3% [23]. The IRIS curve (blue 
line) shows the flux-weighted redshift distribution, dl /dz, of IRIS sources from the CIB model of 
Ref. [37], used for the fitting described in Section 3.2. Each distribution has been normalised so that 
the area under it is unity. 




[36]. The catalogue consists of 13,823 clusters of galaxies of richness 10 < N < 190 over 
redshifts 0.1 < z < 0.3 in an area about 7500 deg 2 , and is estimated to be more than 90% pure 
and more than 90% complete for clusters more massive than 2 x 1O 14 /i" 1 M . 2 Fig. 1 shows 
the redshift distribution of the clusters, compared to the IRIS source distribution used in the 
model described in Section 3.2. See Ref. [23] for further details on the maxBCG catalogue. 

In our work, we construct a maxBCG overdensity map by creating a blank map of zeros 
at the same resolution as the IRIS map (NSIDE = 1024), and then setting to unity any pixel 
containing a cluster coordinate from the maxBCG catalogue. The mean is then removed and 
the map scaled by the inverse number of illuminated pixels, i.e.: 



, m — (m) 

m'= \ 2.1 

where m is the map of ones and zeros and (m) is the mean value of m after applying the 
maxBCG mask (see below). 

Each of the IRIS and maxBCG maps are masked when computing the cross-spectrum 
(Sec. 2.2). For maxBCG, we simply mask all regions falling outside the maxBCG survey. 
The IRIS mask excludes (a) regions with high Galactic infrared emission and (b) resolved 
point sources. For (a), we do a conservative cut and mask all IRIS tiles which have a mean 
map value more than 2.0 MJysr" 1 . This is to minimise potential noise due to Galactic cirrus 
emission. Additionally, we mask a finger-shaped region near the north pole where the IRIS 
map is of particularly low quality. For component (b) , we mask all entries in the IRAS Point 
Source Catalog [34] and the IRAS Small-Scale Structure Catalog [38] with discs of radius 5'. 
The objects in these IRAS catalogues are typically bright (> l Jy), very nearby sources that 
we do not expect to add to our cross-correlation signal. The sky coverage after both IRIS 



2 iV ga i 00 is the number of galaxies within a radius containing a galaxy density which is 200 times larger 
than the average density. 
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Figure 2. The masked, full sky IRIS map, centred on the Galactic north pole with Mollweide's 
projection. The figure shows both the masks — the IRIS mask, to remove high Galactic emission and 
point sources, and the maxBCG mask, selecting the survey area (which is mainly in northern Galactic 
latitudes). The finger-like portion of the mask near the north pole is to remove a flaw in the IRIS 
map. Total coverage is 6.41%, or 2644 dcg 2 . 

and maxBCG masks are applied is 6.41%, or 2644 deg 2 , and contains 4784 clusters. Figure 2 
shows the IRIS 100 |om map with each of the masks. 

2.2 Computation of the Cross-Spectrum 

The cross-spectrum and autospectra are calculated using PolSpice [39]. We correct for a 
Gaussian beam with FWHM of 4.3' [22] as well as for a residual transfer function that is left 
over by PolSpice due to imperfect correction for the 5' point source mask. This transfer 
function, which is small compared to the uncertainty in our power spectrum measurement, 
was determined by measuring the power spectrum of a synthetic map of white noise after 
applying the same mask. Finally, we bin the cross-spectrum into bands of width A£ = 150 
below I = 2000; at higher multipoles, where the data are noisier, we choose a bin width of 
A£ = 450. The cross-spectrum is shown in Fig. 3. 

To estimate our uncertainties, we generate 144 random realisations of maps with the 
same power spectrum as the IRIS map, and then cross-correlate this map with the cluster 
map in the same way as described above. Their average spectrum, consisting of the mean of 
these 144 synthetic spectra in each A£-band, is consistent with a null- measurement, having 
X 2 = 18.4 for 13 degrees of freedom (dof), with a probability to exceed (PTE) of 0.14. This 
null-test is shown in Fig. 3. We take the covariance of these mean band powers as our 
estimate of the measured spectrum uncertainties. 

The choice to mask IRIS panels with mean intensity greater than 2.0MJysr~ 1 is con- 
servative, as is the 5'-radius masking of IRAS point and compact sources. Nevertheless, in 
order to verify that the masks do not bias our result (which could potentially occur if, for 
example, fewer maxBCG clusters are detected where the Galactic emission is brightest), we 
computed power spectra with different mask sizes. We varied both the mask size based on 
the IRIS mean intensity, as well as the radius of point-source masks. Table 1 lists all the mask 
combinations that we performed along with their respective sky-coverage. Figure 4 shows 
that the different mask sizes do not significantly alter the cross-spectrum. In the end, the 
2.0MJysr~ 1 IRIS threshold plus 5' point source mask combination was chosen as it balanced 
sky area against low cross-spectrum uncertainties. 
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Figure 3. IRISxmaxBCG cross-spectrum. The cross-spectrum bandpowers are shown in the upper 
panel along with a power-law fit with Ct oc ^~ 128 , as well as the best fit that uses the phenomenological 
IRIS DSFG model from Ref. [37], with a parametrisation of the cluster bias (see Section 3.2). The 
lower panel shows a null test taken by cross-correlating random-realisations of maps having the IRIS 
autospectrum with the maxBCG map; note that this plot has a linear scale on the y-axis. The result 
is consistent with zero power. 



3 Model fitting 

3.1 Power law fit 

We first model the IRISxmaxBCG cross-spectrum as a power law of the form Cp = A (£/1000) 7 
This simple model fits the data well, with x 2 pe r dof = 8.2/13 (PTE = 0.83). We measure 
A = (1.3 ± 0.2) x 10~ 7 MJy sr^ 1 and 7 = -1.28 ± 0.12, where the uncertainty in A includes, 
and is dominated by, the 13.5% uncertainty in the absolute flux calibration of the DIRBE 
maps used to re-calibrate the original IRAS maps by Ref. [22] . 

A model assuming no correlation between the IRIS and maxBCG overdensity maps, 
where the cross-spectrum is zero at all scales, yields a \ 2 °f 221. Comparing to the power- 
law fit, we infer that IRISxmaxBCG cross-correlation is detected at a significance level of 
over 40cr. 

3.2 Large-scale cluster bias 

We also use our cross-spectrum measurement to constrain the large-scale bias of the maxBCG 
clusters, in order to compare with previous clustering analyses using the maxBCG catalogue. 
To do this, we use an existing model for the properties of the unresolved dusty galaxies in 
the IRIS map. Ref. [37] (hereafter A12) present a phenomenological model for the luminosity 
function, spectral energy distribution and clustering properties of CIB sources, constrained 
by differential number counts and angular power spectra from 70 \vm. to 1.4 mm. We write 
the IRISxmaxBCG cross-spectrum at multipole I as (following, e.g., Refs. [11, 40]): 

Ci = J (1 \ z) (W*' z ))( fe ciust(fc, z)) (j„(z))cut j ^7^~7^7 PDM ( fc ' *)L-'+V 2 » ( 3 - 1 ) 



- 5 - 



Point Source 
Mask (arcmin) 



2 
5 
10 



IRIS Intensity Threshold (MJysr ) 
1.65 2.0 2.5 3.0 



6.57% 
6.41% 
6.00% 



10.34% 14.46% 



Table 1. Sky coverage of masks used to test for systematics. The overall mask was varied by 
admitting different thresholds of IRIS average intensity (columns), while the point-source mask was 
varied in radius (rows). The fiducial mask, used for all our analysis, is the central entry in bold type. 



1-3 



CU- 
I 

a* 



Varying IRIS mean intensity mask 



cut 



Varying point source mask radius 



Fiducial 2.0 MJysr 



le-06 



le-07 



le-08 



Fiducial 5' Radius 
2' Radius 
10' Radius 




1000 



1000 



Multipole (I) 



Figure 4. Tests for systematic errors due to the choice of masks (c.f., Table 1), which show that 
our mask choice is not biasing the cross-spectrum. Left: The mean intensity cut for allowing IRIS 
tiles into the unmasked region is varied from 1.65 to 3.0MJysr _1 . Right: The radius of point-source 
masks is varied from 2 to 10'. 



where % is the comoving distance to redshift z, b ga \ is the mean IRIS source bias, 6 c i us t is 
the mean cluster bias, j u is the IRIS source emissivity density, dN c \ ust /dz is the redshift- 
distribution of the maxBCG clusters and Pdm is the linear matter power spectrum. Since 
our spectra are from scales £ > 100, we apply the small-scale Limber approximation [41, 42], 
setting k = (£+ 1/2) j\- The subscript 'cut' denotes the fact that bright sources (Siooum ^ 
1 Jy) have been masked in the IRIS map; expressions for (6 ga i) and {j u ) C ut are given in 
Section 2 of A12. Fig. 1 shows a comparison of the redshift distributions of the model's IRIS 
sources and the maxBCG galaxy clusters. 

We repeated the fit described in A12 with the addition of the IRISxmaxBCG cross- 
spectrum. We found that the cross-spectrum cannot be described using a scale- independent, 
linear, cluster bias; this is not surprising since over the range of physical scales probed by 
our measurement (0.15 < fc/Mpc -1 < 7.5), the matter power spectrum becomes strongly 
non-linear. Part, or all, of the non-linear clustering may arise from IRIS sources that are 
actually members of the maxBCG clusters (the 'one-halo' term of the halo model clustering 
formalism; for a review see Ref. [43]). In the present analysis, we do not attempt to separately 
model this term, and instead simply parametrize the cluster bias relative to the linear matter 
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power spectrum using a simple phenomenological expansion in powers of k: 



(bcinst(k,z)) = b 



( y q ) 



(3.2) 



with two free parameters, the large-scale bias, bo, and a parameter describing the scale 
dependence, ^4bias- The pivot wavenumber ko is chosen to be 1 Mpc" 1 . We find that the 
data are unable to constrain higher-order terms in this expansion. 

Compared to the original A12 fit, the x 2 increases by Ax 2 = 11.2 for 13 extra degrees 
of freedom (15 IRISxmaxBCG bandpowers minus two new parameters, bo and ^bias)- This 
corresponds to a PTE of 0.59, indicating that the IRISxmaxBCG cross-spectrum is well- 
described by the A12 CIB source model with the addition of the cluster bias parameters. 
None of the A12 model parameters are significantly shifted compared to the fit without the 
IRISxmaxBCG cross-spectrum. The large-scale maxBCG cluster bias is constrained to be 
bo = 2.6 db 1.4, which agrees with the measurement of 2.80 ± 0.13 from Ref. [44], obtained 
using the real-space maxBCG cluster correlation function. The nuisance parameter ^4bias 
is measured to be 3.9 ± 2.5, and is highly anti-correlated with bo (see below). The best-fit 
IRISxmaxBCG model is plotted in the top panel of Figure 3. 



4 Discussion 

While we have measured a strong cross-correlation signal, our constraint on the cluster bias 
is weak compared to the real-space maxBCG correlation function measurement of Ref. [44]. 
This is primarily because, as stated, our measurement probes non-linear scales, and so is 
unable to distinguish the large-scale cluster bias from non-linear effects, leading to a strong 
degeneracy between bo and ^4bias- Extracting the quasi-linear and linear-scale correlation 
(k < 0.1 Mpc -1 , corresponding to £ < 100 for this redshift range) using IRIS data is diffi- 
cult because of effective noise from Galactic dust emission. We expect the imminent release 
of submm and microwave Planck maps to lead to significant improvements in this mea- 
surement. While the Planck resolution is comparable to IRIS, having multiple frequency 
channels will both facilitate reducing noise from Galactic dust and provide a stronger test of 
CIB source models. A joint analysis of the large-scale cluster correlation function and our 
cross-spectrum measurement, with an appropriate treatment of the covariance, would also 
help to disambiguate the linear and non-linear clustering. 

It is intriguing that the shape of the IRISxmaxBCG cross-spectrum is similar to that of 
the clustered dusty galaxy component measured from existing FIR, submm and microwave 
measurements (Cg oc £~ -L2 , e.g., Ref. [16, 31, 45]). Future work is clearly required to 
understand the implications of the apparent similarity, and the extent to which the small- 
scale CIB-cluster cross-correlations can be used to tell us about the occupation of clusters 
with dusty galaxies or how star- forming regions are distributed within clusters. The Her- 
schel SPIRE instrument has channels at 250, 350 and 500 (a.m and roughly ten times better 
resolution than IRIS or Planck, and should provide improved constraints on CIB-cluster 
cross-correlations on smaller angular scales. SPIRE sky coverage includes over a hundred 
square degrees overlap with SDSS and upcoming optical surveys. 3 

The significant correlation between the IRIS map and the maxBCG clusters qualitatively 
supports the prediction of a moderate degree of spatial correlation between the thermal 

' s http: //www . astro . caltech. edu/~ viero/viero_homepage/helms_hers .html 
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Sunyaev Zel'dovich effect associated with groups and clusters and faint dusty galaxies that 
make up the CIB [46] , relevant for current efforts to constrain the kinetic Sunyaev Zel'dovich 
effect [31]. We note, however, that the maxBCG clusters are all at z < 0.3, whereas most 
of the tSZ power comes from groups and clusters at higher redshift (see, e.g., Figure 1 of 
Ref. [47]). Nevertheless, our result suggests that useful quantitative constraints on the tSZ- 
CIB correlation may be obtained by cross-correlation of catalogues of galaxy clusters at 
higher redshifts with FIR, submm and microwave maps. 

We have not accounted for any spatial offset between BCG position and cluster centre in 
our analysis. Ref. [48] find that such offsets mostly lie in the range —2.5 < log(Ar//i _1 Mpc) < 
—1.0, where Ar is the physical offset of the BCG from the cluster centre of mass, using strong 
gravitational lensing. While such offsets would lead to a suppression of small-scale power in 
the IRISxmaxBCG cross-spectrum, their effect is small for £ < 3000. We similarly have 
not accounted for impurity in the maxBCG catalogue. Clearly, these uncertainties may be 
expected to become important for analyses attempting a more physical interpretation of 
CIB-cluster cross-correlations. 

5 Conclusions 

We have measured the cross-power spectrum of the 100 (xm IRIS map and an overdensity 
map constructed from the maxBCG catalogue of optically-selected clusters at 0.1 < z < 
0.3. The cross-spectrum is detected at a statistical significance of over 40o\ This result is 
possible because the dominant components of the IRIS auto-correlation power spectrum — 
Galactic dust on large scales and Poisson-like power from self-correlation of DSFGs on small 
scales — are not present in the cross-spectrum signal, allowing information about the spatial 
correlation of the faint IRIS sources and maxBCG clusters to be cleanly extracted. 

While our result provides only a fairly weak test of the A12 model of CIB sources, it 
has demonstrated the potential of CIB-optical cross-correlations. We anticipate rapid im- 
provements in this field with upcoming analyses using Planck and Herschel data. Looking 
beyond our analysis, which has focussed on galaxy clusters, a wealth of information may 
be gleaned using other optically-selected populations (e.g., quasars and luminous red galax- 
ies). CIB-Quasar cross-correlations, for example, may help elucidate the connection between 
galactic nuclei activity and star formation, and how this relationship evolves with redshift. 
The use of large catalogues of objects with well- measured redshifts in cross-correlations will 
also provide important constraints on the redshift-distribution of unresolved CIB sources. 
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